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Abstract 

Theories which include neutrinos in large extra dimensions can be constrained 
by nuclear beta decay experiments. We examine universality of (5 decay strengths of 
Fermi transitions. From this we find that the extra dimensional Yukawa coupling for 
a higher dimensional scale of 10 TeV, and two extra dimensions can be constrained to 
y < 1. Kinematic implications are also discussed. In particular, an extra dimensional 
scenario will produce a tritium decay endpoint spectrum with a different shape than 
that for just one massive state. 
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1 Introduction 



It is well known that consistent string theory can only be formulated in space time higher 
than four dimensions, specifically 10 or 11 dimensions. Until recently the predominant 
view was that since the extra dimensions have not been detected, they must be compact- 
ified in small volumes with radii of Planck size. This was supported by the study of the 
heterotic string where the string scale Ms is related to the Planck scale Mp via |l j 

M s = Mpsffo^ (1) 

where «g denotes the unified gauge coupling constant and k is an integer of order unity 
which labels the level of the Kac-Moody algebra. Here the string scale is close to the 
Planck scale and this renders string effects unobservable. Recently, Witten || observed 
that in nonperturbative string theory Ms can be much lower. Lykken |5| and Antoniadis 
[f§ independently suggested the possibility that Ms can be as low as a few TeV. In this 
case, the fundamental string scale is not constrained theoretically as above but instead is 
constrained phenomenologically. This gains support from breakthroughs in understanding 
D-brane constructions in string theory ||. In the D-brane scenario the number of extra 
dimensions scanned by the Standard Model (SM) particles and the graviton can be very 
different. 

Buoyed by these developments intense efforts are made to understand the main features 
of low scale string theories. A popular construction is to assume a factorizable geometry. 
The gravitational field is taken to propagate in the full 10 or 11- dimensional bulk volume 
whereas the SM particles are localized in the 3-brane and hence are not sensitive to the 
extra dimensions 0, ||, || 01- The crucial observation |J is that the weakness of gravity 
in four dimension is due to its spreading in the n extra dimensions and Eq.(l) is replaced 
by 

Ml = M: +2 V n (2) 

where V n is the volume of the compactified extra space and M* is the higher dimensional 
Planck scale. The two scales Ms and M* are related but the exact relation will not be 
needed in this paper. It is simplest to consider toroidal compactification, for which the 
volume is 

V n = {2 7 r) n R 1 R 2 R 3 ...Rn (3) 

and Ri(i — 1, 2, . . . n) are the radii of the extra dimensions. For symmetric compactifica- 
tion one takes all the radii to be equal and denoted it by R. 

This formulation has dramatic effects that lead to a new perspective on the mass 
hierarchy problem || and the gauge coupling constants unification ||. Furthermore, 
these theories predict a deviation from the Newtonian 1/r law of gravitational potential 
at submillimeter range. The fact that this law is well tested above a distance of 1mm 
leads to the limit that n > 2 for = 1 TeV, for the symmetric case. New experiments 
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currently underway will probe smaller distance scales ||. Besides the direct search for 
a deviation of the Newtonian gravitational law the above scenario can also be indirectly 
constrained by astrophysical considerations. Although somewhat model dependent these 



considerations suggest that n > 2 and M* > 30 TeV JlO . 

A different scenario to explain the weak gravitational force is proposed by Randall 
and Sundrum [IT|. This mechanism does not require large extra dimensions but instead 
invokes a non-factorizable geometry. Here the extra dimension has the effect of modifying 
the space-time metric. Experiments that probe gravity at small distances are then not 
interpreted in terms of the scaling law of Eq. 0. The phenomenology of the graviton and 
its Kaluza-Klein (KK) excitations is different from the case of factorizable geometry and 
has been studied in Ref.(T^|. Given the differences of the these two scenarios and their 
extensions, and the richness of the phenomenology for each case it is clear that probes 
other than the graviton and its KK spectrum will be important additions to the study of 
higher dimensional physics. 



It has been suggested recently [O, that an extra right-handed neutrino that is 



a standard model gauge singlet may also be a probe of the extra dimensions. When a 
field is allowed to extend into the compactified extra dimensions it has associated with 
it Kaluza-Klein excitations which can have detectable effects in low energy experiments. 
Neutrinos propagating in the extra dimension will be referred to as the bulk neutrinos. 
The couplings of the bulk neutrino states with the active neutrino will give rise to a small 
neutrino mass for the active neutrino due to the spreading in the bulk volume. Some 
effects of these neutrinos in various low energy observables have been suggested A 
solution to the solar neutrino problem in terms of matter enhanced flavor transformation 



of v e L into bulk neutrino states was studied in [16|. Bulk neutrinos can induce energy 



loss in stars and the possible connection to the supernova collapse phase is studied 
in 
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In this paper we concentrate on the constraints on the extra dimensional scenarios from 
universality tests of nuclear /3-decays. Such studies have been used to provide a strong 
basis for the SM and more recently their usage is mainly in the test of unitarity of the 
quark mixings. We point out that the difference in Q-values for the nuclei which are very 
well measured can be used to constrain the properties of the bulk neutrinos such as the 
Yukawa couplings, the higher dimensional scale, M*, and the compactification radius R. 
These are very general properties of neutrinos in extra dimensions and one does not need 
to employ the scaling relation between the Planck scale and M* as given in Eq.(fJ). To 
illustrate the effects of the bulk neutrinos we give the Kurie plots of two nuclei and the 
nuclear recoil spectrum of a Fermi transition. The simplest model of bulk neutrinos given 
below may allow a signal to be detected in the 3 H beta decay spectra. Other models which 
modify the tritium beta decay endpoint are discussed in However, in searches for a 
small admixture of a more massive neutrino, the features of the extra dimensional model 
are not within currently detectable limits. However more realistic models with larger 
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mixings with the active neutrinos may be constructed which can be probed by these 
experiments as well. We concentrate only on one family of bulk neutrino and the active 
u e L- However, it is straightforward to extend the model to include three bulk neutrinos; 
one for each family. This implies further model assumptions and will take us into the 
realm of models for fermion families in extra dimensions which is beyond the scope of this 
study. 

This paper is organized as follow. In section ^| we briefly outline the model of bulk 
neutrinos that we use and the mixings with v e ^. This serves the purpose of fixing of our 
notation. Sec. ^| details the study of universality using the data from Q- values of several 
/9-decays. Next we present the impact of bulk neutrinos on the beta spectrum, and look 
at the recoil momentum for Fermi transitions. We give our conclusion in Sec [5| 



2 A Simple Model of Bulk Neutrinos 

The basic idea of producing neutrino masses using extra dimensions is given in Ref. 



T6| . In the interest of being self contained and to establish our notation we give some 
of the construction of a simple model of bulk neutrinos. The simplest model assumes 
that the fermions charged under the standard model gauge group as well as the gauge 
bosons and the Higgs boson are localized on a 3-brane embedded in the bulk of larger 
dimensions. If we assume only one extra dimension then space-time is labelled by (x 11 , z) 
where /i = 0, 1, 2, 3. The extra coordinate z is assumed to compactify into a circle of radius 
R. The brane scenario stipulates that the left-handed SM lepton doublet L = {y eL ,tj}) is 
given by L(x M , z = 0). The fields v e L and cl are Weyl spinors. Next we assume that a SM 
singlet fermion, denoted by u(x, z) exists and it propagates in the full five dimensional 
bulk. Its right- and left-handed projections are labelled by vr and respectively. It is 
distinguished from the SM neutrino by not carrying a flavor index. In this simple model 
gravity is weak because there are additional dimensions in which it propagates. 

In five dimensions there are five gamma matrices r M and T 5 where /i = 0,1,2,3. A 
convenient representation of the Clifford algebra in five dimensions is to choose T M = 7 M 
and T 5 = ry 5 with the usual 7^ of Minkowski space-time. The effective Lagrangian for 
generating a neutrino mass is given by 

r2irR r2nR 

C= dzv {i^da + iT 5 d z ) v + y* dz5(z)LHu R + h.c. (4) 
Jo Jo 

where is the dimensionful Yukawa coupling and H is the Higgs doublet and it is related 
to the dimensionless coupling y via 

and n = 1 for one extra dimension. 
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For simplicity we have neglected a higher dimensional bare Dirac mass term. This 
can be naturally implemented under Zi orbifold compactification [14|j. Following the 
Kaluza-Klein ansatz we Fourier expand v R and vi as follows: 

Substituting this into Eq.(|]) and integrating over z yields the effective Lagrangian in four 
dimensions: 

oo oo oo 

£= PkLil^d^ kL + v kR i^d ll v kR }+ m kVkL"kR+ —/== LHv kR + h.c. (8) 

k=~ oo fc=— oo V ^VTix k= _ 00 

where ^ 

mk = (9) 

is the mass of the k — th KK tower state. The mass splitting between each adjacent 
tower state is 1/R. As seen in Eq. (§) the coupling between the KK-tower neutrino states 
and the active v e ^ given by the Yukawa term. After spontaneous electroweak symmetry 
breaking a Dirac mass term is generated for the active v e ^ and is given by 

m D = - yV (10) 
where v = 247 GeV. This mass is suppressed by a bulk volume factor as first noticed in 



Ref. 131 and 14 



Rewriting all the mass terms in Eq. ( |8|) we have 

oo oo 

m D v eL v 0R + m D Y^ v eL (i> kR + v_ kR ) + fn k {v kL v kR - V- kL V-kR) + h.c. (11) 

k=l k=l 

We find it useful to define the following orthogonal states 



\r = -j= Ofcfl + V-hr) "In = -^m(y kR - u_ kR ) (12) 



and 

v 'kL = ~^ UkL ~ U - kL ^ V kL = ~^ VkL + V - kL ^ ( 13 ) 

which can be used to cast Eq.ll into the form 

oo oo 

m D u eL u 0R + V2m D ^ v eL v kR + ^ m k {u kL u' kR + v' kL v kR ) + h.c. (14) 

k=l k=l 
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The states with double prime superscripts have no low energy interactions and will be 
ignored. The mass terms can now be written in the familiar form of vlMvr in the bases 
of u L = (u eL , v' 1L , v' 2L , . . .) and u R 
states looks as follows 



v 'qri v '\ri v 2Ri • • •)• The mass matrix M for k + 1 



M 



( m D \[2m D \[2m D 

o h 





2_ 

R 



y/2m D \ 





R / 



(15) 



\ o o oo 

To find the left mass eigenstates, we consider the matrix MM^ . Explicitly, 



/ (k + \)C C 2C ... K \ 



1 



V 



c 

K 



l 
o 
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Ok 2 ) 



(16) 



where ( = \f2moR. The mass eigenvalues A are given by the characteristic equation of 
Eq.(|l^) which is det [MM^ — A 2 ]. After some algebra it reduces to the transcendental 



equation |L4 



7rC 2 cot(7rA J R) = 2XR 



(17) 



The solution of this equation yields the lowest mixed mass eigenstate ~ mo- For small ( 
this is also the maximally mixed state. It is reasonable to assume that experiments that 
directly measures the neutrino mass will be probing this state. The current experimental 
limit on the mass of v e is m Ve < 2.5 eV p0[ . When combined with this, Eq. (|I~0"D gives 
the following constraint 



y<3.2x 10 



-3 



AC. 



lOTeV 



R \ 2 fm D \ 
1mm/ VleV/ 



The diagonalization of Eq. 
eigenstate Uq L given by ||1^ 



]) leads to the mixing of the fc-th state to the lowest mass 



tan2#k 



k 2 -{k + \)C 



(19) 



This equation diagonalizes the submatrix between the lowest state v§ L and the state v kL . 
Where ( is small, this is a good approximation to the exact mixing angle, obtainable by 
complete diagonalization of the matrix MMl In this limit, ( can be used as an expansion 
parameter in the diagonalization process. For example, for = 10 TeV, y = 10~ 3 and 
R ~ 10~ 7 mm we find £ ~ 10 _1 . Hence, for small ( we can see from above that the 
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mixing of the KK states with the lowest mass neutrino eigenstate becomes progressively 
smaller as the value of k increases. Furthermore, these KK states have no direct gauge 
interactions with other SM particles and their presence can only be probed through their 
mixing with the state v™ L . 

In the small £ limit, the vacuum oscillation survival probability for electron neutrinos 
is a particularly simple expression, 



k 2 



1 - 2 sin" 



2E„ 



+ 0(( 4 



(20) 



where 8m 2 , 



Ok 



m D — m 



k- 



The data from a reactor neutrino experiment can be used to 



find a limit on the Yukawa coupling. CHOOZ [BTJ quotes a limit on two neutrino mixing 
of sin 2 26 < 10% for 5m 2 > 10~ 3 eV 2 . This translates into a limit on the extra dimensional 
parameters of 

o/lmm\ 1/2 / \5 
, < 4.4 x 10- (_) (— ) . (21) 

If combined with the limit from the neutrino mass, Eq [L8|, we find that the least stringent 
limit on y occurs at R = 7 x 10~ 6 mm. 

However, Eq.(|l9"D reveals a different phenomenon if ( is large; i.e. ( ^> 1. Now it is 
the higher KK states with k « £ that will have the largest mixing with the the lowest 
eigenstate; whereas states with lower k values will have small mixing of order ^. In order 
for this to happen in n = 1 case we require R > 10~ 3 mm. Obviously this cannot happen 
with £ < 1 and a small compactification radius. When relatively large values of R are 
considered, masses of these KK states can be in the ~ 1 eV range. 

The above considerations can be generalized to higher extra dimensions. For n = 2 
and different compactification radii Ri and R2 the scaling equation for mjj becomes 



yv 



2nM*^2Rjr 2 



[22] 



and yields the following constraint the on the Yukawa coupling 

( M» W R X R 2 



y<1.8x 10 b 



VlOTeV/ Vlmm 2 



2 ( mD \ 

viev; 



(23) 



For Ri = R 2 and a comparison with Eq.flTBD shows that this is a much less stringent 
constraint on the parameters of the theory compared to the case of only one extra di- 
mension. This is expected since the bulk neutrino state has a larger volume in which 
to spread and the Yukawa coupling is inversely proportional to the square root of this 
volume. An interesting case occurs for asymmetric compactification where one radius is 
much smaller than the other. Take for example R 2 = 10 _11 mm and R\ = 0.1mm then 
Eq.(^) approaches the case of n — 1 and gives a limit on y similar to that of Eq.(|l8|). 
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The masses of the KK tower states are generalized to 



\ 




(24) 



where k and / denote the KK level corresponding to radii Ri and i? 2 respectively. It is 
seen that the smaller compactification radius R2 will give rise to higher mass KK state 
in the asymmetric compactification scenario. For example if R2 = 10~ 7 mm one can have 
keV bulk neutrinos mixing into the active v e state that come from the small radius (i.e. 
for k = and small I values). When both radii are large such neutrinos come only from 
the high KK values. Next we examine the mixings of these bulk states. 



For small mixing the electron neutrino state can written as |L6 

•>. = & £ r« + <■ s jmwr^ ) (25) 

where Ci = y/2m£ ) R 1 « 1 and N is the normalization constant. It is seen that for 
the asymmetric case the KK states which come solely from the small radius may have 
larger masses than those from the larger radius; moreover, they are accompanied by 
correspondingly smaller mixings. We conclude that high mass states will have small 
mixing, regardless of the size of the compactification radii mainly due to the constraint 
onrnc from the direct limit on the electron neutrino mass. Since the direct limits on the 
mu and tau neutrino masses are much less strict, large mixing may occur in these families 
for considerably smaller KK masses. 

For two extra dimensions, the constant, N, depends on the cutoff scale M*, as iV ~ 
1 + ( 2 log(M*i2), if Ri = R'2 — R. Limits derived from universality of Fermi transitions 
which are presented in the next section do not depend on this normalization and therefore 
on the cutoff procedure. However, other calculations such as vacuum oscillation survival 
probabilities will depend on the normalization. In particular, this would become impor- 
tant if one were to compare the theory with the reactor neutrino data. In this case, the 
constraint on y and M* from universality could be used together with survival probability 
data to limit the size, R. 

In the next section we study the constraint (3 decay universality places on the param- 
eters y, M* and R. 



3 Universality 

In this section we study the consequences for beta decay of the model for bulk neutrinos 
outlined in section |[ Although bulk neutrinos do not have gauged interactions with 
SM particles, the allowed Yukawa coupling enables them to generate mass for the active 
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neutrino. As seen in the last section the smallness of the neutrino mass is due to the 
spreading of the state in the higher dimensions. The KK excitations are also seen to have 
mixing effects with the active neutrino. In a given weak decay, the number of KK states 
that can be excited will depend on the energy released in the reaction. If 1/R « Q, 
where Q is the nuclear Q-value, then many states contribute. For example, for n = 1, if 
Q = lMeV and R = mm then around 500 states contribute. In previous studies of 
such classical weak interactions one concentrates on one or two massive neutrino states 
and their kinematic effect. In the present scenario, towers of KK neutrinos ranging in 
mass from eV to a few MeV are involved in a nuclear (3 decay. They have to be taken into 
account in the study of universality, electron spectra and nuclear recoil spectra. We shall 
examine these issues separately to see if they provide useful constraints on the Yukawa 
coupling, y and the size R, of the extra dimension. We begin with universality. 

The rate of beta decay from a single state in the parent nucleus to another state in the 
daughter nucleus is given by, 

A = |PK) (26) 

The ft- value is proportional to the inverse of the matrix element. P(m v ) is the phase 
space factor, which depends on the mass of the emitted neutrino, m u . We define the phase 
space factor as 

P(m v ) = — / F(Z, E e )E ePe (Q + m e - E e )[(Q + m e - E e f - ml} 1/2 dE e (27) 

The Q-value is the energy difference between the initial and final nuclear states, while E e , 
and p e are the electron energy and momentum respectively. The Coulomb wave correction 



factor is denoted by F(Z, E). We use the approximation of this factor found in [22 |. 

For pure Fermi, + — > + transitions and zero mass electron neutrinos, the corrected 
ft- values, Ft = /t(l + 5^)(l — 5c) for all nuclei should be the same under the assumption 
of universality in the standard model. Here, 5r is the nucleus dependent part of the 
radiative correction and 5c is the isospin symmetry breaking correction. These corrections 
contribute at the 1% level. Results of Ft- values for 10 nuclei are listed in ||23|| . 

The bulk neutrino scenario outlined in the previous section would cause an apparent 
deviation from the constant Ft values predicted by universality in the standard model. 
This is because the actual phase space factor in an extra dimensional scenario is more 
complicated than Eq.([27|). In general ft- values are found by using Eq. [26|, with a measured 
rate and a calculated phase space factor. The transition rate for the case of one dimension 
is given in the small mixing limit, ( << 1, by 



In 2 

A 



xd 




p( mi/ «0)+£p(m !/fc «-)| 5 



(2* 
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The maximum number of neutrinos to be summed over is determined by the maximum 
mass of a neutrino that can be released in the beta decay, k max = QR. Transitions with a 
higher Q value will have more neutrino states that can contribute to the decay. Note that 
the ft value in an extra dimensional scenario will not necessarily take on the same value 
as given in the standard model. (Although, the 'true' Ft- values will be the same for all 
+ — > + transitions.) This will change the calculated value of the quark mixing angle V ua - 
as discussed below. Because of this uncertainty, it is best to employ a normalization when 
looking for the effects of bulk neutrinos. In this example we will compare the relative 
difference in beta decays from two different nuclei. 

Assuming that bulk neutrinos exist in one extra dimension, we can compare the ft- 
value for the extra dimensional scenario with the apparent /i-values for the standard 
model for the nucleus, A±, by taking a ratio, 

fJfAU „ L *K 2 \ k f? P(A u m Uk ) c 2 (2Q) 

SM, apparent 

where we have used Xsu/^xd = 1- Since in the extra dimensional scenario, the ft- values 
for two nuclei are the same up to the isospin and radiative corrections, the apparent 
Standard Model values will not be the same, 

ft(A 2 )sM,apparent _ , fc ™^ R P(A U m Uk ) ( 2 _ k ^=^ R P (A 2 , m„J ( 2 



SM, apparent fc tt P(A,0) k 2 ^ P(A 2 ,0) k 2 

-MAj) - 6 R (A 2 )] + [5 C (A 1 ) - 5 C {A 2 )\. (30) 

The maximum value of the second two terms may be obtained from the experimental 



uncertainty in the Ft- values from | 23]| . The sums in the above expression are very insen- 
sitive to lower bound on the sum, so it is only necessary to consider the highest modes in 
the sum, where there are relatively large differences in the phase space factor ratio. The 
mixing of the low energy KK states is not relevant to universality considerations. 

The above formula is valid for large (, provided that ( 2 /k 2 << 1, if the term (( 2 /k 2 ) is 
replaced by f (() 2 {( 2 / k 2 ) . Here f(() is a function which must be calculated for each £ by 
explicit diagonalization of the mass matrix given by Eq. (|l6|) . This explicit diagonalization 
shows that the mixing of state with k >> m^R is approximately ('/k = f{()(/k. For 
( = 100, /(C) = 0.007. For larger £ /(() decreases, while for C « 1, /(C) = 1- 

We use f3 + decay data from 14 and 54 Co with a combined error of 0.22%. These nuclei 
have Q values of 1.81 MeV and 7.22 MeV respectively. This can be translated into a limit 
on the Yukawa coupling,?/ and the size of the extra dimension, R, 

10TeV\ /lmm\ 



0.0022 > 1.9 x 10- 7 MeVV/(C) 2 (- 



R ) 



Q ' R /p(A u m Uk ) P(A 2 ,m Uk )\ 1 <W P(A 2 ,m Vk ) 1 



t\ P(A U 0) P(A 2 ,0) )ml k k ^ iR P(A 2 ,0) ml k 



(31) 
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Converting the sums in the above equation into integrals and integrating over all the 
available neutrino states produces a limit on the Yukawa coupling. As long as the small 
mixing limit applies, the limit on y has little R dependence and remains fairly constant 
for example, at y < for a scale of M* = 10 TeV. The limit becomes less strict for 
a larger high dimensional scale. The lack of R dependence can be seen when converting 
the sum in Eq. fl3T|) to an integral, and making the change of variables dk —>■ Rdm v . In 
the case of large mixing, the limit is less strict, due to the factor /(C)- 

When comparing with Eqs.(|18|, pT]) , one can see that the limit on the Yukawa coupling 



from universality is always less stringent than the limit on the Yukawa coupling derived 
from the maximum neutrino mass and mixing determined by tritium beta decay and 
reactor neutrino experiments, respectively. This latter limit is shown as the bottom curve 
in Figure 1. Therefore, one does not obtain a useful constraint from universality when 
considering only one extra dimension. 

One can perform the same exercise for two extra dimensions which both have the same 
size R. In the perturbative limit, the mass eigenstates are mj tn « (k 2 + n 2 )/R 2 . In this 
case the sum ranges over two indices, k and n. An equation similar to Eq.([H]) may be 
derived for this case, 



« 9 9 /10TeV\ 2 /mm\ : 
0.0022>6xlO-MeW(— ) (— ) 



/ f(J,,m, t J _ P(A,,m„ k J \ 1 _ P(A 2 ,m 

^ \ T>( A. <W P( A~ C\\ J ™2 2^ 



M=i v P (^'°) P (^,0) k ^ QlR P(A 2 ,0) m^j 



(32) 



All two dimensional cases considered here, M* > ITeV, fulfill the condition C 2 << 1- 
We solve the above equation by converting the sums to integrals, and find that the limit 
is again insensitive to R. Figure 1 also shows the constraint on M* and y for the two 
dimensional case. It can be seen that greater precision is needed to probe scales greater 
than 10 TeV, since above this scale the Yukawa coupling is not constrained to be less than 
of order 1. 

The Standard Model value of V u d can be obtained from measured Ft values as in ]23[ 
with the additional input of Gf, the weak coupling constant from muon decay. Using 
these quantities, the unitarity sum V 2 d + V 2 S + V 2 b is calculated to be smaller than one 
by two standard deviations. The extra dimensional scenario would effect both the muon 
decay measurement as well as the Ft values from beta decay. This introduces additional 
parameters and model uncertainties, such as the coupling of the muon neutrino to the 
same or another Kaluza Klein tower of neutrinos. If, for example, the muon neutrinos did 
not couple to any bulk neutrinos, then the calculated value of V u d should appear too large 
and an extra dimensional scenario could not account for the shortfall in the unitarity sum 
through the determination of V u & from beta decay. 
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Next we make contact with another much discussed test of universality using charged 
pion decays f24 |. The pion decays, 7r — > ev e and tt — > /iz/^ can be probed kinematically. 



Since the electron in the decay is relativistic, while the muon is not, the contributions 
of a massive neutrino differ considerably in the two decays; a massive electon neutrino 
will make a larger difference than a massive muon neutrino. In the case of massless 
neutrinos the standard model ratio is T(7r — > ez/)/T(7r — > fiu) = 1.233 x 10 -4 . The 
experimentally measured quantity is = (1.230 ± 0.004) x 10~ 4 Therefore, the 
extradimensional contribution to these decays can not exceed 0.3%. The contribution 
has been calculated explicitly by |24| and for the n = 2 case with two families of bulk 



neutrinos. The contribution is approximately given by 

, » ZlTeV\ 2 , ZlTeV\ 2 9 ^ 

7.527T x 10Ye (^-J - 6.237T x 10- 1 (— j < 0.003. (33) 

With this reaction alone, the only way to obtain a constraint on either y e or y^ as a function 
of M* is to make an assumption about one of the couplings such as, y e ~ 1 or ~ 1. 
Ideally, one would obtain a separate constraint on one of the couplings and use the pion 
decay data to limit the other. Taking the constraint on y e from beta decay universality 
discussed previously we find that another two orders of magnitude in precision in the 
beta decay data would be required in order to probe y^. However, a complete discussion 
must include an analysis of neutrino oscillations which is beyond the scope of this paper; 
however, see p6| for a discussion. 



4 Kinematic searches 



In addition to changing the apparent ft values, bulk neutrino scenarios can in princi- 
ple also be probed kinematically. Here we demonstrate the consequences of the extra 
dimensional scenario on nuclear recoil momenta, and electron spectra. Massive neutrino 
searches in beta decay have focused on examining beta spectra |20|, |27|, ^8], ^9|] and on 
examining the nuclear recoil spectra [ |3"0"| , |3T|. A tower of Kaluza-Klein bulk states will 
produce a somewhat different signal than just one massive neutrino. We examine the 
signature in both of these types of spectra from the KK tower of states. Larger mixing 
scenarios may be probed by tritium beta decay searches for eV range neutrinos. However, 
as we shall see below, the model of extra dimensions cannot currently be constrained by 
searches for heavy (keV-MeV) neutrinos. 

With a massive electron neutrino, the number of electron neutrinos as a function of 
momentum would be, 



N( Pe 



p 2 e F(Z,p e 



1 1 

ml ft 



(Q - E e + m e y -ml (Q - E e + m e ) 



1/2 



(34) 
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We show the normalized Kurie plot, (N{p e )/{p 2 ,F(Z, Pe))) 1 ^ 2 vs. E e for 38m K at the top 
of Figure |3], for the case of m u = 0. The Q- value for this decay is 5.02 MeV. 
In the case of one extra dimension, the function is modified to 



N( Pe 



p 2 e F(Z,p e ) ml ft 



xd k=0 



U, 



fee 



1 1/2 

{Q-E e + m e f -ml (Q-E e + m e ). (35) 



Here |£4 e | 2 is the mixing of the kth mass eigenstate with the electron neutrino. For heavy 
massive neutrino searches (keV-MeV) very small mixing is relevant as discussed below. 
However, for massive neutrino searches in the eV range, such as tritium beta decay, a 
slightly larger mixing is applicable. Several states will give a different signature than just 
one or two massive states in a kinematic study. We illustrate this in the lower panel of 
Figure 0. This panel shows the beta endpoint spectrum for tritium for three cases. The 
dashed line shows the spectrum for a single massive neutrino of 2.3 eV. The solid line 
shows the result for a scenario with one extra dimensions, and £ = .13, 1/R = 25 eV giving 
mo = 2.3 eV. The first ten mass eigenstates were calculated by explicit diagonalization 



of the matrix, Eq. (16). The first three occur at 2.3 eV, 25 eV and 50 eV, with mixings 
\Uke\ 2 of 0.974, 0.017 and 0.004 respectively. As can be seen from the figure, the Kaluza- 
Klein tower produces several bumps as well as less counts slightly away from the endpoint 
than a single massive neutrino. For comparison a two neutrino mixing scenario is shown 
as the dot-dashed line. The mass eigenstates are at 2.3 eV and 25 eV with a mixing of 
sin 2 6 « 2.5%. 

We now turn to searches for massive neutrinos in the keV and MeV range. In the case 
of one extra dimension and small mixing, the function is modified to 



N{p e 



1 1 



P 2 e F(Z,p e ) ml ft. 



xd 



1--^) {Q-E + m e ) + 



kmax /" 2 

J2 j^\(Q-E e + m e ) 2 -m 2 
k=i fc 



1/2 



(Q - E e + m e 



(36) 



We have normalized to the number of counts in order to eliminate the uncertainty in ft x d- 
Therefore, we replace N(p e ) — > N(p e )/N, where for one extra dimension, N = f N(p e )dp e . 
We plot the results for the extra dimensional scenario as the ratio, 



n 



N{p e 



1 1/2 



Np 2 F(Z,p) 



I 



- xd 



N( Pe 



11/2 



Np 2 F(Z,p) 



(37) 



J SM 



Two cases are shown in Figure The panel in the bottom left hand corner shows a 
situation where many modes can contribute to the decay; R = 10 _8 mm, m^ = leV. The 
mass of the lightest mode is at about 0.02 MeV and has a mixing of around ( 2 = 5 x 10 -9 . 
Nuclei with smaller Q values will have a slightly larger signature for the continuum case. 
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The experimental limits on the mixing angle for a single 17 keV neutrino, for example 
27| are of order 10~ 3 . Since the mixing of the ~ 20 keV modes from this extra dimensional 



scenario are considerably weaker, these modes are not likely to be detectable with present 
experimental data. The mixing of a ~1 MeV neutrino is smaller still, since the mixing 
goes as (/k = \/2m D R/k. If the lowest KK state occurs around 1 MeV, this implies a 
size R ~ 2 x 10~ 10 mm and a mixing of ( 2 ~ 3 x 10~ 13 . The point here is that the mixing 
of the a high mass KK mode is proportional to which is constrained by tritium beta 
decay experiments. Therefore, in this simple model the mixing of keV and MeV neutrinos 
will always be small. The bottom right hand corner of Figure |3] shows contributions from 
three separate neutrino masses, 1, 2 and 3 MeV. The bumps shown in the lower right 
panel are unobservable, since current constraints on the mixing of an MeV neutrino are 
many orders of magnitude less (sin 2 29 ~ 10~ 3 ). 

We also consider the method of detecting massive neutrinos by measuring nuclear 
recoil spectra. In the presence of a massive neutrino, the nuclear recoil spectrum for a 
pure Fermi transition has the shape, 

j r Emax 

P r (Z,r) = - F(Z,E e 



a ' 2 2 2^ 



rE e E u + r-(r - p e - p 



dE e (38) 



Here, r is the recoil momentum of the nucleus and a = 1. For a pure Gamow- Teller 
transition, a = —1/3. The limits of integration depend on the mass of the neutrino and 
are given by 



e =i 

J - / max,mm ^ 



E (Eq — r 2 + me 2 — m 2 ) ± rJ (E 2 — r 2 — m 2 — ml) 2 — km 2 v m 2 e 



El - r 2 



(39) 



Here Eq = Q + m e is the total energy available for the decay. We again normalize the 
distribution by integrating over all possible recoil energies and finding P r (m u ), for each 
given neutrino mass. The distribution for m v = is shown on the top panel of Figure 
The normalized probability distribution for the scenario with one extra dimension looks 
like, 

pj, = (1 - ^f-)P r (r,m u « 0) + Ek £Pr(r,m u « k/R) 

P r 

where the normalization is 

P r = (1 - ^p) J P r {r, m v w 0)dr + £ £ J P r {r, m v w k/R)dr (41) 

k 

To see the effect of the bulk neutrino states on recoil spectra, we again plot normal- 
ized ratios, (P r (r)/ P r )xd/{P r {r) / P r )sM The results are shown in Figures |] for the Fermi 
transition in 38m K. The top panels show the normalized spectra for a single zero mass neu- 
trino. The bottom left panels show the effect when many neutrino modes can contribute 
to the decay, while the bottom right panels show the case where only a few neutrinos 
can contribute. However, as in the case with the Kurie plots, the mixing of the KK 
modes is several orders of magnitude too small to be detected with current experimental 
measurements. 
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5 Conclusions 



We have studied the effects of bulk neutrinos in nuclear (3 decays in a simple model which 
only allows the SM singlet state to propagate in one or more compactified extra dimen- 
sions. All other SM particles are confined to a 3-brane and have no KK excitations. There 
are many well studied nuclei with different Q values that allow us to use universality to 
place constraints on the Yukawa coupling and the fundamental scale M*. This informa- 
tion is complementary to the various studies of physics of higher dimensions using the 
gravitational force or KK graviton probes. 

Our study shows that for n — 1 the most stringent constraint on the parameters of 
the bulk neutrinos arise from the direct measurement of the mass of v e and from reactor 
neutrino oscillation experiments. Universality tests do not give additional information. 
On the other hand for n = 2 the universality test restricts the Yukawa coupling to be 
less than unity for M* < 10 TeV. Future experiments with radioactive beams measuring 
decays with Q values of order 10 MeV and improved accuracy with current measurements 
that can push this by a factor of ten will be most welcome. 

The Kurie plot and the recoil nuclear spectrum in some selected decays are also given. 
For probes of massive Kaluza-Klein states we find the nuclei with smaller Q values are 
more sensitive. Currently, for medium and heavy nuclei, the sensitivity of experiments 
designed for massive neutrino searches are several orders of magnitude below what is 
needed to probe mixing of the massive KK neutrinos in these reactions. This may be 
due to the simplicity of the model we have considered and hopefully such studies will 
prove useful for more realistic models with enhanced mixings that also address the family 
problem. We also point out that our study may be used in connection with other data to 
constrain parameters of second and third generation mixing with bulk neutrinos. These 
scenarios are more model dependent and contain more parameters. We give an illustration 
of this point for charged pion decays. 

On the other hand we find that the tritium (3 decay spectrum can be significantly 
altered by the presence of KK bulk neutrino states in the eV range. They can reveal 
themselves in the shape of endpoint region of the Kurie plot if the mixings are favorable. 

It is useful to compare the results from universality limits to the ones derived from 
astrophysical constraints which are indirect. The latter relies on the energy loss carried 
away by bulk neutrinos. For example with limits derived from the neutrino magnetic 
moment |L7] for n = 2 there is essentially no constraint on the neutrino magnetic moment 
induced from the sort of extra dimensional scenario presented here, while the limit on the 
Yukawa coupling is of order 1 for M* ~ 10 TeV from universality. For n = 1 one obtains a 
constraint from the mu and tau neutrino magnetic moment whereas universality gives no 
useful constraints. As mentioned before in this case for electron neutrinos the strongest 
constraint comes from the electron neutrino mass limit and oscillation experiments. 
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Figure 1: The upper curve shows the limit on the extra dimensional Yukawa coupling 
from nuclear beta decay as a function of higher dimensional scale M*. This constraint is 
for a number of extra dimensions, n=2. The lower curve is for n=l. It shows the limit on 
the extra dimensional Yukawa coupling from the tritium beta decay neutrino mass limit 
and the CHOOZ oscillation probability limit. This lower curve is a tighter constraint than 
the universality test for one extra dimension. 



Figure 2: The top panel shows a the endpoint of a Kurie plot for tritium beta decay, 
for which the Q value is 18.6 keV. The horizontal axis is electron kinetic energy - Q. The 
bottom panel shows the same region with three different neutrino mass scenarios. This 
plot shows the subtraction of zero mass neutrino Kurie plots from the massive neutrino 
Kurie plots. The dashed line shows the spectrum for a single mass neutrino of 2.3 eV. 
The solid line shows the spectrum for one extra dimension with ( = .13 and 1/R — 25 eV. 
The dot-dashed line shows a two neutrino mixing scenario with states at 2.3 eV and 25 
eV and a mixing of 2.5% 



Figure 3: The top panel shows a normalized Kurie plot for 38m K. The horizontal axis 
is total electron energy. The dashed line corresponds to a zero mass neutrino, while the 
solid line corresponds to a neutrino with mass 2.3 eV. The bottom two panels show ratios 
of the normalized Kurie plot in the standard model to the normalized Kurie plot in the 
extra dimensional scenario. The bottom left panel, shows the effect for R = 10 _8 mm 
and a Dirac mass term of m D = 1 eV. The bottom right panel was produced with the 
parameters R = 2 x 1CT 10 mm and m D = 1 eV. 



Figure 4: The top panel shows a normalized nuclear recoil spectrum for 38m K. The 
bottom two panels show ratios of the normalized spectrum in the standard model to the 
normalized spectrum in the extra dimensional scenario. The bottom left panel, shows the 
effect for R = lCT 8 mm and a Dirac mass term of m D = 1 eV. The bottom right panel 
was produced with the parameters R = 2 x 1CT 10 mm and mn = 1 eV. 
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